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Abstract 



I Charged Current reactions of Ue on ^^C have been studied using a fi'^ 

decay-at-rest Ue beam at the Los Alamos Neutron Science Center. The 
cross section for the exclusive reaction ^^C(i'e, e~)^^Afg.3. was measured 
to be (8.9 ± 0.3 ± 0.9) x 10~^^ cm^. The observed energy dependence of 
the cross section and angular distribution of the outgoing electron agree 
well with theoretical expectations. Measurements are also presented for 
inclusive transitions to ^'^A'^ excited states, ^^C{ue, e~y^N* and compared 
with theoretical expectations. The measured cross section, (4.3 ± 0.4 ± 
0.6) X 10"*^ cm^, is somewhat lower than previous measurements and 
than a continuum random phase approximation calculation. It is in better 
agreement with a recent shell model calculation. 



*now at University of Alabama, Tuscaloosa, AL 35487 
^now at Indiana University, Bloomington, IN 47405 
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1 Introduction 



In recent years neutrino interactions with nuclear targets have been used to 
detect low energy neutrinos (< 500 MeV) from many sources: solar, atmo- 
spheric, supernova explosions, reactors and accelerators. An understanding of 
the nuclear cross sections is necessary for interpretation of the measurements. 
Especially important nuclei are ^^C and ^^O because of the widespread use of 
oil and water detectors. 

Many calculation techniques have been used to determine neutrino-nuclear 
cross sections. Shell model techniques work best at lower energies where transi- 
tions to continuum states are not large. At intermediate energies the Continuum 
Random Phase Approximation (CRPA) is frequently used, while at still higher 
energies the Fermi gas model is thought to work well. Comparison of different 
calculations of the cross section for a particular process can provide an indi- 
cation of the uncertainty involved. Experimental measurements of some cross 
sections are, however, necessary to establish the range of validity of the different 
calculation techniques. 

Relatively few measurements of neutrino-nucleus cross sections exist in the 
energy region of the present experiment, Ei, < 52.8 MeV. The best mea- 
sured nucleus is carbon for which three experiments, including the Liquid Scin- 
tillator Neutrino Detector (LSND) experiment, have previously reported re- 
sults. E225Q at LAMPF, the KARMEN Collaboration § at the ISIS facil- 
ity of the Rutherford Laboratory and LSND[^ have measured the cross sec- 
tion for the exclusive reaction ^^C{ve^e''Y'^ Ng,s, and for the inclusive reaction 
-^^C(i/e, e~)^^iV* to all other accessible ^"^N final states. The ^"^N ground state 
reaction dominates the total yield as it is the only allowed {I = 0) transition 
that occurs in this process. The cross section for producing the ^"^N ground 
state can be calculated to an accuracy of approximately 5% as it can be rep- 
resented in terms of form factors jj] that can be reliably extracted from other 
measurements. Calculation of the inclusive cross section for transitions to ex- 
cited states of ^^A^ is much less straightforward. Various theoretical techniques, 
each with their own strengths and limitations, have been used to calculate the 
cross section 1^, ^ 0, |[ ^ ^ 0. Comparison with measurements may help 
clarify the theoretical picture. In this paper we report our final results for these 
processes, including measurements of the angular distribution of the electron 
with respect to the u^, direction and the energy dependence of the ground state 
transition. 

Measurements also exist for two processes closely related to carbon scat- 
tering: ^~ capture on ^^Cjl2[ and scattering on carbon using a beam of 

from 7r+ decay-in-flight (DIF)||l^. Because these three processes occur at 
different energies, i?^, and momentum transfers, Q, they constrain different as- 
pects of theoretical calculations. A good test of a theoretical procedure is its 
ability to predict all three processes. For the Ve carbon measurement Ei, « 32 
MeV, Q « 50 MeV/c and the inclusive cross section is dominated by transitions 
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to low multipoles (1^,1 ,2 ). In contrast, for the v^j_ carbon measurement 

« 180 McV, Q « 200 MeV/c and excitatfons occur up to 100 MeV. The 
capture process, which occurs from the S state, is intermediate between these 
two processes with Q « 90 MeV/c. 

The measurement of the inchisive cross section for ^"^Civ^, ^~Y'^N* 
several years ago by LSND attracted substantial interest because a CRPA 
calculation[H predicted a cross section almost twice as large as that observed. 
An improved calculation by the same group 0| together with an improved calcu- 
lation of the neutrino energy spectrum and flux has reduced but not eliminated 
the discrepancy with the measured cross section. Recent calculations using the 
shell model § | are in better agreement with the measured cross section. 

Hayes and Towner|P calculated a cross section of 4.1 x 10"''^ cm^ for the 
process ^^C(i/e, e~)^^A^* using the same shell model procedure that provided 
the best agreement with data for the reaction ^^C(z^^, /i~)^^A^*. This is lower 
than both the earlier 1^ and the more recent CRPA calculations by Kolbe et 
al. of 6.3 X lO""*^ cm^ and 5.5 x 10"'*^ cm^ respectively. Thus measurements of 
this cross section can provide a useful test of the relative merits of the different 
theoretical techniques that have been used. 

2 The Neutrino Source 

The data reported here were obtained between 1994 and 1997 at the Los Alamos 
Neutron Science Center (LANSCE) using neutrinos produced at the A6 proton 
beam stop. We chose to exclude data obtained in 1998 from this analysis because 
only electrons with reconstructed energies above 20.4 MeV were fully processed 
for that year's data. Since this analysis is dominated by systematics we decided 
to limit our analysis to the 1994-1997 data which has a uniform efficiency above 
18 MeV. The neutrino source is described in detail elsewhere [p^. In 1994 and 
1995 the beam stop consisted of a 30 cm water target surrounded by steel 
shielding and followed by a copper beam dump. The high-intensity 798 MeV 
proton beam from the linear accelerator generated a large pion flux from the 
water target. The flux of v^, used for the measurements reported here arise from 
the decay at rest (DAR) of stopped 7r+ and /x"*". This decay chain yields almost 
equal intensities of Ve, V/i and Vfj, with the well-determined energy spectra shown 
in Fig. I. 

After the 1995 run the beam stop was substantially modified for accelerator 
production of tritium (APT) tests. The most significant change for the analysis 
presented in this paper was the replacement of the water target by tungsten 
and other materials with high atomic number. This resulted in reduced 7r+ 
production and a lower DAR neutrino flux, largely due to the change in the 
neutron to proton ratio in the target. 

The corresponding decay chain for tt^ and fi~ is highly suppressed due to 
three factors. First, production of vr^ is smaller than for 7r+. Second, tt^ which 



3 



10 20 30 40 50 
Neutrino Energy (MeV) 



Figure 1: Flux shape of neutrinos from pion and muon decay at rest. 

stop are absorbed by nuclear interactions. Finally, most /i^ which stop are 
absorbed before they can decay. These stopped /i^ arise from tt^ which decay 
in flight. 

A few percent of the 7r+ (tt^) produced at the beam dump decay in flight 
to Vfj, (vfi) with energies up to 300 MeV. Those t'^ above muon production 
threshold provide the source used for our measurement of ^^C(f^, ^~)^^A^*]l3|. 
In the analysis of the reaction ^^C{h'e, e^y^N* presented in section ^ below we 
correct for a small background arising from /i^ produced by high energy and 

The LANSCE beam dump has been used as the neutrino source for previ- 
ous experiments [p^ p^ , pTf . A calibration experiment p8[| measured the rate of 
stopped /i^ from a low-intensity proton beam incident on an instrumented beam 
stop. The rate of stopped /i+ per incident proton was measured as a function 
of several variables and used to flne-tune a beam dump simulation program ||l9| . 
The simulation program can then be used to calculate the flux for any par- 
ticular beam dump configuration. The calibration experiment determined the 
DAR flux to ±7% for the proton energies and beam stop configurations used at 
LANSCE. This uncertainty provides the largest source of systematic error for 
the cross sections presented here. The LANSCE proton beam typically had a 
current of 800 /iA at the beam stop during the 1994-1995 running period and 
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Figure 2: Detector enclosure and target area configuration, elevation view. 

1000 fiA for 1996-1997. For 1994 and 1995 the energy was approximately 770 
MeV at the beam stop due to energy loss in upstream targets, while it was 
approximately 800 MeV in 1996 and 1997. The water target was out for 32% 
of the 1995 data. Upstream targets contributed 1.4% to the DAR flux in 1994 
and 1995. The DAR i^e flux averaged over the LSND detector was 3.08 x 10" 
cm-2 for 1994 and 3.45 x 10" cm-^ for 1995. 

The 1996-1998 data were obtained with the APT beam stop. There were 
no upstream targets for almost all of the data taking for this period. The DAR 
i^e flux averaged over the LSND detector was 1.32 x 10^"^ cm^^ for 1996 and 
2.73 X 10"'^^ cm~^ for 1997. For the full data sample used in this paper the i/^ 
flux is 10.58 X 10" cm-2. 

3 The LSND Detector 

The detector is located 29.8 m downstream of the proton beam stop at an 
angle of 12° to the proton beam. Figure || shows a side-view of the setup. Ap- 
proximately 2000 g/cm^ of shielding above the detector attenuates the hadronic 
component of cosmic rays to a negligible level. The detector is also well shielded 
from the beam stop so that beam associated neutrons are attenuated to a neg- 
ligible level. Enclosing the detector, except on the bottom, is a highly efflcient 
liquid scintillator veto shield which is essential to reduce contributions from the 
cosmic ray muon background to a low level. Reference provides a detailed 
description of the detector, veto, and data acquisition system which we briefly 
review here. 

The detector is an approximately cylindrical tank containing 167 metric 
tons of liquid scintillator and viewed by 1220 uniformly spaced 8" Hamamatsu 
photomultiplier tubes (PMT) covering 25% of the surface inside the tank wall. 
When the deposited energy in the tank exceeds a threshold of approximately 

4 MeV electron-equivalent energy and there are fewer than 4 PMT hits in the 
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veto shield, the digitized time and pulse height of each of these PMTs (and of 
each of the 292 veto shield PMTs) are recorded. A veto, imposed for 15.2 ^s 
following the firing of > 5 veto PMTs, substantially reduces (10"'^) the large 
number of background events arising from the decay of cosmic ray muons that 
stop in the detector. Activity in the detector or veto shield during the 51.2 /xs 
preceding a primary trigger is also recorded, provided there are > 17 detector 
PMT hits or > 5 veto PMT hits. This activity information is used in the 
analysis to reject events arising from muon decay. Data after the primary event 
are recorded for 1 ms with a threshold of 21 PMTs (approximately 0.7 MeV 
electron-equivalent energy). This low threshold is used for the detection of 2.2 
MeV 7's from neutron capture on free protons. The processes measured in this 
paper, ^^C(i'e, e~Y'^Ng,s. and ^^C(z/e, e~)^^iV*, do not produce neutrons. Thus, 
in the present analysis, detection of 2.2 MeV 7's is used to help determine beam- 
related backgrounds with associated neutrons. The detector operates without 
reference to the beam spill, but the state of the beam is recorded with the event. 
Approximately 94% of the data is taken between beam spills. This allows an 
accurate measurement and subtraction of cosmic ray background surviving the 
event selection criteria. 

Most triggers due to electrons from stopped muon decays (Michel electrons) 
are identified by a preceding activity produced by the decay muon. Occasionally, 
however the muon will not satisfy the activity threshold of > 17 detector PMT 
hits or > 5 veto PMT hits. For several LSND analyses, including the present 
one, it is desirable to further reduce the number of unidentified Michel electrons. 
Therefore, for data recorded after 1994 all PMT information was recorded for 
a period of 6 /xs (2.7 muon lifetimes) preceding certain primary events. For 
the present analysis this "lookback" information is used to further reduce the 
cosmic ray muon background as described in Section |^. 

The detector scintillator consists of mineral oil {CH2 ) in which is dissolved 
a small concentration (0.031 g/1) of b-PBD[^. This mixture allows the sep- 
aration of Cerenkov light and scintillation light and produces about 33 photo- 
electrons per MeV of electron energy deposited in the oil. The combination of 
the two sources of light provides direction information for relativistic particles 
and makes particle identification (PID) possible. Note that the oil consists al- 
most entirely of carbon and hydrogen. Isotopically the carbon is 1.1% ^^C and 
98.9% 12c. 

The veto shield encloses the detector on all sides except the bottom. Addi- 
tional counters were placed below the veto shield after the 1993 run to reduce 
cosmic ray background entering through the bottom support structure. More 
counters were added after the 1995 run. The main veto shield|^ consists of a 
15-cm layer of liquid scintillator in an external tank and 15 cm of lead shot in 
an internal tank. This combination of active and passive shielding tags cosmic 
ray muons that stop in the lead shot. A veto inefficiency < 10~^ is achieved 
with this detector for incident charged particles. 
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4 Analysis Techniques 



Each event is reconstructed using the hit time and pulse height of all hit PMTs 
in the detector jl^. The present analysis relies on the reconstructed energy, 
position, direction, and two PID parameters, xj^^ and a, as described later 
in this section. The particle direction is determined from the Cercnkov cone. 
The parameters x'tot ^'^'^ ^'^s used to distinguish electron events from events 
arising from interactions of cosmic ray neutrons in the detector. We directly 
measure the response of the detector to electrons and neutrons in the energy 
range of interest for this analysis by using copious control data samples. We 
also use a GEANT Monte Carlo simulation, LSNDMC(2^ to simulate events in 
the detector. 

The response of the detector to electrons was determined from a large, es- 
sentially pure sample of electrons (and positrons) from the decay of stopped 
cosmic ray /i* in the detector. The known energy spectra for electrons from 
muon decay was used to determine the absolute energy calibration, including 
its small variation over the volume of the detector. The energy resolution was 
determined from the shape of the electron energy spectrum and was found to 
be 6.6% at the 52.8 MeV end-point. 

There are no tracking devices in the LSND detector. Thus, event positions 
must be determined solely from the PMT information. The reconstruction pro- 
cess determines an event position by minimizing a function Xr which is based on 
the time of each PMT hit corrected for the travel time of light from the assumed 
event position to the PMTjl^. The procedure used in several previous analyses 
systematically shifted event positions away from the center of the detector and 
thus effectively reduced the fiducial volume pSf. The reconstruction procedure 
has been analyzed in detail and an improved reconstruction procedure was de- 
veloped which reduces this systematic shift and provides substantially better 
position resolution. This procedure also provides results which agree well with 
positions obtained from the event likelihood procedure described in Ref. |Q. 
In the analysis presented in this paper, a fiducial cut is imposed by requiring 
13 > 35 cm, where D is the distance between the reconstructed event position 
and the surface tangent to the faces of the PMTs. Events near the bottom of 
the detector (y < —120 cm) are also removed, as discussed in Section ^. 

The particle identification procedure is designed to separate particles with 
velocities well above Cerenkov threshold from particles below Cerenkov thresh- 
old. The procedure makes use of the four parameters defined in Ref. |p^ . 
Briefly, Xr ^nd Xa are the quantities minimized for the determination of the 
event position and direction, xt is the fraction of PMT hits that occur more 
than 12 ns after the fitted event time and Xtot is proportional to the product of 
Xr, Xa and Xt- 

Several previous LSND analyses |l^, |2^ have used Xtot for particle iden- 
tification. The distribution of Xtot for electrons, however, has a small variation 
with electron energy and with the position of the event. Therefore, in this pa- 
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Figure 3: Particle identification parameters (a) xj^j and (b) Xa for electrons 
and neutrons. In the present analysis we require Xa < 4.0 as indicated by the 
arrow in (b). 



per, we used a modified variable, x^^j, with a mean of zero and sigma of one, 
independent of the electron energy and positions. We also used the variable, 
a, which is based on the event likelihood procedures discussed in Ref. 
As in ]2^ , a comes from a separate event reconstruction than that which pro- 
duced x'tof It is similar to the parameter p discussed there, which is based 
on the ratio of Cerenkov to scintillator light. The a parameter varies from 
to 1 and peaks at one for electrons and at 0.3 for neutrons. The combination 
Xa — Xtot + 10(1 ^ Q^) provides better separation of electrons and neutrons than 
Xtot by itself. 

Figure ||(a) shows the x'tot distribution for electrons from stopping /i decay 
and for cosmic ray neutrons with electron equivalent energies in the 18 < i?e < 
50 MeV range. Neutrons, after thermalizing, undergo a capture reaction, n + 
p ^ d + J. The 2.2 MeV 7's are used to select a clean sample of neutron 
events. For a neutron E^, is the equivalent electron energy corresponding to the 
observed total charge. Figure ||(b) shows the Xa distribution for the same events. 
In the present analysis we eliminate most cosmic ray neutron background by 
requiring Xa < 4.0. We note that a modest particle identification requirement 
was imposed in the initial data processing that created the samples analyzed 
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here. The effect of this requirement is also included in the analysis. 

Beam-off data taken between beam spills play a crucial role in the analysis 
of this experiment. Most event selection criteria are designed to reduce the 
cosmic ray background while retaining high acceptance for the neutrino process 
of interest. Cosmic ray background which remains after all selection criteria have 
been applied is well measured with the beam-off data and subtracted using the 
duty ratio, the ratio of beam-on time to beam-off time. The subtraction was 
performed separately for each year's data using the measured duty ratio for 
that year. The ratio averaged over the full data sample was 0.0632. Beam-on 
and beam-off data have been compared to determine if there are any differences 
other than those arising from neutrino interactions. Any differences are small 
and the 1.1% uncertainty in the duty ratio introduces a negligible effect in the 
present analysis. 

5 Electron Selection Criteria 

In this section we describe the selection criteria used to obtain a clean sample 
of inclusive electrons arising from neutrino interactions in the detector. In the 
next section we present the analysis of the relatively pure sample of events from 
the process ^^C(z/e, e~Y'^Ng,s. which we obtain by requiring the detection of the 
positron from the /3-decay of the ^'^Ng,s.- Section |^ then presents the analysis 
of the reaction ^'^C{ve, e~Y'^N* using the sample of inclusive events without an 
identified positron from the /3-decay of ^^Ng,s. 

A lower limit of 18.0 MeV is imposed on the electron energy to eliminate 
the large cosmic ray background from ^^B /?-decay as well as most 15.1 MeV 
gamma rays from the neutral current excitation of carbon. The ^^B nuclei arise 
from the absorption of stopped on ^^C nuclei in the detector. The scattered 
electron from the reaction ^"^Ciy^, Ng g has a maximum kinetic energy of 
35.5 MeV due to the Q value of 17.3 MeV. Allowing for energy resolution we 
impose an upper limit of 40 MeV on the electron energy. 

The selection criteria and corresponding efficiencies for electrons with 18 
MeV< Ee < 40 MeV are shown in Table |l]. The reconstructed electron position 
is required to be a distance D > 35 cm from the surface tangent to the faces 
of the PMTs. The requirement y > —120 cm removes a small region at the 
bottom of the detector for which the cosmic ray background is relatively high 
due to the absence of a veto below the detector. There are 3.34 x lO'^" ^^C 
nuclei within this fiducial volume. The fiducial volume efficiency, defined to be 
the ratio of the number of events reconstructed within the fiducial volume to 
the actual number within this volume, was determined to be 0.880 ±0.055. This 
efficiency is less than one because there is a systematic shift of reconstructed 
event positions away from the center of the detector as discussed in Section ^ 

Several selection criteria are designed to further reject cosmic ray induced 
events. Events with more than three veto PMT hits or any bottom counter 
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Table 1: The electron selection criteria and corresponding efBciencies for events 
with 18 MeV< < 40 MeV. 



Quantity 


Criteria 


Efficiency 


Fiducial volume 


D > 35 cm, 


0.880±0.055 




y > —120 cm 




Particle ID 


Xa < 4 


0.940±0.018 


In-time veto 


< 4 PMTs 


G.988±0.010 


Past activity 


See text 


G.635±0.012 


Future activity 


Atf > 8.8 ^s 


0.991±G.003 


Lookback 


likelihood 


0.994±0.004 


DAQ and tape dead time 




0.962±0.010 


Total 




0.492±0.035 



coincidence during the 500 ns event window are eliminated. The past activity 
cut is designed to reject most electron events arising from cosmic ray muons 
which stop in the detector and decay. This background has a time dependence 
given by the 2.2 /is muon lifetime. The past activity selection criteria reject 
all events with activity within the past 20 fis with > 5 veto PMT hits or > 17 
detector PMT hits. We also reject any event with a past activity within 51 
US with > 5 veto PMT hits and > 500 detector PMT hits. A small (0.5%) 
loss of efficiency arises from a cut (made during initial data processing) on past 
activities that are spatially correlated with the primary event, within 30 /zs of 
the primary event and have > 4 veto PMT hits. 

Muons which are misidentified as electrons are removed by requiring that 
there be no future activity consistent with a Michel electron. Any event with 
a future activity with fewer than 4 veto PMT hits and more than 50 detector 
PMT hits within 8.8 /iS is rejected. 

Cosmic ray muons which fire < 6 veto PMTs (10^"^ probability) and stop 
in the iron walls of the detector will not register as past activities. Some of 
the decay electrons will radiate photons which will enter the detector and be 
reconstructed as electrons within the fiducial volume. In previous analyses we 
simply relied on the beam-off subtraction procedure to remove this background. 
Here we use the "lookback" information described in Section || to reject events 
from this source. This results in slightly smaller statistical errors in the final 
beam excess sample. 

For primary events with > 300 PMT hits and no activity within the past 
35^s (20^s) for 1995 data (1996-1998 data), we recorded all PMT information 
for the 6 fis interval proceeding the event. Muons with < 6 veto PMT hits 
will appear in this "lookback" interval as a cluster of veto PMT hits spatially 
correlated with the primary event. The distribution of time between the veto 
signals and the primary event should be consistent with the muon lifetime, and 
the distributions of veto PMT hits and veto pulse height should be consistent 
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Figure 4: The distribution of time between the primary and the veto signal 
for beam-off events rejected using the "lookback" information compared with a 
curve corresponding to the muon hfetime. 

with that measured for muons producing < 6 veto PMT hits. We developed 
a likelihood procedure based on these distributions which allowed us to reduce 
the beam off background by 9% with only a 0.6% loss of efficiency for neutrino 
events [p5| . Figure |^ shows the time between the veto signal and the primary 
for rejected events. The fitted lifetime of 2.08±0.07 fis agrees well with the 
expected average lifetime of 2.12 /is for stopping /x"*" and /i~ in oil. 

The acceptances for the past activity, the future activity, the "lookback" 
and the in-time veto cuts are obtained by applying these cuts to a large sample 
of random events triggered with the laser used for detector calibration. These 
laser events are spread uniformly through the run and thus average over the 
small variation in run conditions. The acceptance for the 15.1 /is trigger veto is 
included in the past activity efficiency. 

A sample of Michel electrons was analyzed to obtain the acceptance of elec- 
trons for the PID cut. Figure H] compares the Xa distribution of the inclusive 
electron sample with a Michel electron sample. The agreement is excellent. 
To eliminate any energy dependence, the Michel electrons are given weights 
as a function of energy so that the weighted spectrum agrees with the energy 
spectrum of electrons from the reaction ^'^C{v^,e~'Y^ Ng ^ . The acceptance. 
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Figure 5: The Xa distribution of the beam excess inclusive electron sample. The 
histogram shows the Xa distribution of Michel electrons. 



however, is very insensitive to the assumed energy spectrum. The beam ex- 
cess distribution shown in Figure ^ is obtained by subtracting the beam off 
distribution from the beam on distribution as discussed in Section ||. 



6 The Transition to the ^^TV Ground State 



The reaction + ^^C + ^'^Ng,s. is identified by the detection of the e~, 

followed within 60 ms by the positron from the /3-decay of the ^^iVg.s.. Transi- 
tions to excited states of ^^N decay by prompt proton emission and thus do not 
feed down to the ^^A'' ground state or contribute to the delayed coincidence rate. 
The scattered electron has a maximum kinetic energy of 35.5 MeV due to the Q 
value of 17.33 MeV. The /3-decay has a mean lifetime of 15.9 ms and maximum 
positron kinetic energy of 16.33 MeV. The cross section to the ^"^N ground state 
has been calculated by several groups. The form factors required to calculate 
the cross section are well known from a variety of previous measurements. This 
cross section and the known I'e flux are used to obtain the expected electron 
kinetic energy spectrum. Figure ^ shows the observed and expected electron 
energy distributions for events with an identified /3-decay. Figure shows the 
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Figure 6: The observed and expected (solid line) energy distributions for elec- 
trons from ^^C(fe, e~)^^A^g.s.. 



13 




Figure 7: The spatial distribution of the electron for beam-excess events com- 
pared with expectation (solid line) from ^'^C{iye, e~y^ Ng,s.. 



observed and expected spatial distributions of the same events. Both figures 
show good agreement with expectations. 

Table ^ gives the selection criteria and efficiencies for the ^'^N /3-decay 
positron. Figure ^ shows the observed /3-decay time distribution. The best 
fit curve yields a lifetime of 16.3 ± 0.8 ms in good agreement with the expected 
value of 15.9 ms. Figure ^ shows the distance between the reconstructed elec- 
tron and positron positions for the beam-excess sample. A cut was applied at 
0.7 m, resulting in an acceptance of (99.2 ± 0.8)%. Following an electron pro- 
duced by a neutrino interaction, an uncorrelated particle, such as the positron 
from ^^B /3-decay, will occasionally satisfy all the positron criteria including 
the requirements of time (60 ms) and spatial (0.7 m) correlation with the elec- 
tron. The probability of such an accidental coincidence was measured by using 
the sample of Michel electrons. The inefficiency caused by the 15.1 /xs veto 
and the DAQ dead time are the same as for the electron. Positrons with 4 or 
more in-time veto hits or any bottom veto coincidence are rejected. The energy 
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Table 2: Criteria to select e+ from Ng,s. beta decay and corresponding efficien- 
cies for the reaction ^^C{i'e, e~Y'^N, 



g.s. 



Quantity 



Criteria 



Efficiency 



(5 decay time 
Spatial correlation 
PMT threshold 

Fiducial volume 
Trigger veto 
In-time veto 

DAQ dead time 



52 ^s< t < 60 ms 
Ar < 0.7 m 

> 100 for 1994, 

> 75 after 1994 
£) > cm 

> 15.1 yUS 
< 4 PMTs 



0.974±0.002 
0.992±0.008 
0.856±0.011 

0.986±0.010 
0.760±0.010 
0.988±0.010 
0.977±0.010 



Total 



0.598±0.015 




Figure 8: Distribution of time between the e~ and e+ for beam excess events 
in the ^^C(i^e, e^)^^iVg.s. sample. The best fit curve (solid line) corresponds to 
a lifetime of 16.3 ± 0.8 ms. 
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Figure 9: Distribution of the distance between the reconstructed positions of the 
e~ and e+ for beam excess events in the ^^C(z^e, e~)^^-/Vg.s. sample compared 
with Monte Carlo expectation (solid line). The calculated accidental contribu- 
tion is shown by the dashed line. 
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Figure 10: Observed and expected (solid line) e"*" energy distribution for the 
^"^CiiycC'Y^Ng.s. sample. 

distribution of the positron is calculated from the ^'^N /3-decay using 
dN , ,o 27rr7 

^==Pei^e(^;™..-^;ef X^-^ (1) 

where 77 = Za/ [3e and is the total positron energy (including rest energy) . 
The slight modification of the spectrum due to the shape correction factor [p6| 
was found to have a negligible effect on the results. The ^^iV decays to the 
ground state {E^ax = 16.83 MeV) 94.6% of the time. Beta decay transitions 
to the excited states of carbon are 1.9% {E,nax = 12.39 MeV, followed by a 
4.44 MeV 7), 2.7% {E^ax = 9.17 MeV) and 0.8%{E^ax = 6.5 MeV)[|7). The 
positron annihilates with an electron after stopping. The Monte Carlo was used 
to generated expected distributions for the positron energy and for number of 
hit PMTs. There was a trigger requirement of 100 PMT hits for 1994 data 
and 75 PMT hits after 1994. The beta is required to be less than 18 MeV in 
this analysis. Figure ^ compares the observed and expected positron energy 
distributions. The good agreement shows that the energy calibration is valid 
for these low energy electrons. 

The electron and /3 selection efficiencies, the efficiency for satisfying the elec- 
tron energy requirement and the total efficiency for the reaction ^^C(i^ej e~)^^iVg 
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Table 3: The electron and /3 selection efficiencies, the efficiency for satisfy- 
ing the electron energy requirement and the total efficiency for the process 



Quantity 


Efficiency 


Electron Selection 


0.492 ±0.035 


Future /3 


0.598 ±0.015 


18 MeV < Ee< 40MeV 


0.789 ±0.020 


Total Efficiency 


0.232 ±0.019 



Table 4: Events, efficiency, neutrino flux and flux averaged cross section with 
statistical error only for ^^C(i'e, e^)^^ A^g.s. 



Beam-on events 


743 


Beam-off events 


6 


Accidental background 


4 


Events from ^^C(i^e, e")^^iVg.s. 


733 


Efficiency 


0.232 ±0.019 


i/e flux 


10.58 X 10" cm-2 


{<^) 


(8.9 ±0.3) X lO-'*^ cm^ 



are shown in Table ^ 

Table ^ provides a breakdown of the number of events satisfying the selection 
criteria as well as the total acceptance, the neutrino flux and the resulting flux 
averaged cross section. For the complete data sample the flux averaged cross 
section is (a) = (8.9 ± 0.3 ± 0.9) x lO""*^ cm^ where the ffist error is statistical 
and the second is systematic. The two dominant sources of systematic error are 
the neutrino flux (7%) discussed in Section || and the effective fiducial volume 
(6%) discussed in Section ^. The measured cross section decreases by 1.4% 
when the fiducial volume is reduced by requiring that the electron be at least 
50 cm (instead 35 cm) from the surface of the PMT faces. As discussed in 
section ^ the beam stop was substantially modified after the 1995 run. The 
cross section measured for data taken with the modified beam dump is (5 ± 
6{stat.) )% higher than for the initial beam dump and thus fully consistent 
within the statistical uncertainty. For comparison the previous measurements, 
the final LSND result and several theoretical predictions for the flux averaged 
cross section are presented in Table ^. They are all in agreement with each other. 
We note that all these measurements rely on the neutrino flux calibration that 
is discussed in Section || and thus they have correlated systematic errors. 

For this reaction to the ^'^N ground state it is also straightforward to mea- 
sure the energy dependence of the cross section. The recoil energy of the ^^iV 
nucleus is negligible and thus = + 17.3 MeV where is the electron 
kinetic energy. Figure |ll| shows that the energy dependence of the measured 
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Table 5: Measurements and theoretical predictions of the flux averaged cross 
section for the process ^^C(zyo, e~)^^Ng.s.. 



Experiment 


LSND 


(8.9 ±0.3 ±0.9) X IQ--*^ cm^ 


LSND(previous)|| 


(9.1 ±0.4 ±0.9) X lO-'*^ cm2 


E225(l| 


(10.5 ±1.0 ±1.0) X 10-^2 (,j^2 


karmenH 


(9.1 ±0.5 ±0.8) X lO-'*^ (,ni2 


Theory 


Donnelly^] 


9.4 X lO"''^ cm^ 


Fukugita et al.^ 


9.2 X 10-^2 


Kolbe et al.W\ 


8.9 X 10"*2 


Mintz et al^^ 


8.0 X 10"'*2 



cross section agrees well with expectations [Q. Figure 12 shows the observed 
and expected ^ angular distribution between the electron and the incident 
neutrino. The data agree well with expectations with the jDF = 0.79. The 
mean observed value of cos6' is —0.046 ± 0.021(stat.) ± 0.030(s?/st.) compared 
with an expected value of —0.068. The systematic uncertainty in cos6' is esti- 
mated to be 0.030 based on Monte Carlo studies of the detector response. The 
expected cos Q distribution shown in Figure |l^ includes the effects of experi- 
mental resolution and as a result, is less backward peaked than the theoretical 
distribution used to generate it. 



7 Transitions to excited states of N 

Electrons below 52 MeV are expected to arise from four major neutrino pro- 
cesses: i2c(iyg,e-)i2Ng.,., i2c(i/„,e-)i2N*, 130(1^0, e-)"X and neutrino elec- 
tron elastic scattering. The expected energy and angular distributions of these 
processes are shown in Figures |l3| and |lj, respectively. The different event 
characteristics of these reactions are used to select a sample due primarily to 
the reaction ^'^Ciy^^^e^Y'^'N* . This sample is then used to determine the flux 
averaged cross section and the electron energy and angular distributions for this 
reaction. The measured distributions are compared with expected distributions 
for ^2C(i/e, e~)^^A^* obtained from LSNDMC using the theoretical cross section 
ofRef. 

The reaction ^'^Ciy^^^ e~y^Ng,s. is a source of background since the e+ from 
the /3-decay of ^'^Ng,s. is not always identified. Any event with an identified e"*" 
in delayed coincidence is of course excluded. The background of events with an 
unidentified e"*" is calculated using the positron acceptance given in Table |^ and 
subtracted. 

All three types of DAR neutrinos {i^e, and v^) elastically scatter off elec- 
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Figure 1 1 : The measured and expected (solid line) cross section for the process 

i2C(i.e,e-)i2iV,.,.. 
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Figure 13: The expected distributions of electron energy for (a) ue elastic 
scattering, (b) ^'^C{ve,e-Y'^Ng.s., (c) ^"^Cive^e-y^N* , and (d) ^^C{ve,e-Y^X. 



22 




Figure 14: The expected distributions of cos 9 for (a) ve elastic scattering, (b) 
i2C(!/e,e-)i2Afg.,., (c) ^^C{ue,e-)^^N*, and (d) ^^C{ue,e-y^X. Straight Une 
fits are shown in (b), (c) and (d). 
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trons in the detector but the rate is dominated by scattering. The contri- 
bution due to DIF and scattering on electrons is almost negligible. The 
scattered electron for this process is strongly forward peaked as shown in Fig- 
ure |l^, and thus such events can largely be eliminated with the requirement 
cos6' < 0.9. A measurement of electron-neutrino electron elastic scattering by 
LSND is reported in a separate paper |]30|. 

A third background arises from the interaction of v^. on ^^C nuclei (1.1% of 
the carboiA The expected number of events obtained from the calculated cross 
section[|], ^ for this process is fairly small. The Q value is 2.1 MeV and thus 
about half of the background can be eliminated by requiring an electron energy 
below 34 MeV. We use the cross section calculated by Kubodera|0, 0.525 x 
10~^°cm^, and conservatively assign a 50% uncertainty to this number. The 
KARMEN experiment has measured a cross section of (0.5 ±0.37 ±0.1) x 10^"*° 
cm^ for this reaction[^. 

Most of the excited states of ^"^N decay by prompt proton emission to ^^C 
states and the decays of some of these ^^C states produce gammas between 2 
and 7 MeV. Both these protons and gammas contribute to the visible energy 
detected by LSND and thus increase the apparent electron energy. LSND is rel- 
atively insensitive to low energy protons, partly due to the absence of Cerenkov 
light. For example a 10 MeV proton produces only between 1 and 2 MeV elec- 
tron equivalent energy. The sensitivity to low energy gammas is also somewhat 
lower than for electrons of the same energy. The response to protons, gammas 
and electrons is obtained using events simulated with LSNDMC. For the anal- 
ysis we require that the measured electron energy be between 18 and 34 MeV. 
Including the estimated contributions of protons and gammas based on Ref. 
[ p9| this region contains 72 ±4% of the excited state events. The relatively large 
error assigned is due to uncertainties in the response of LSND, for example to 
low energy protons, as well as estimated uncertainties in the modeling of par- 
ticle production. The upper limit of 34 MeV not only eliminates much of the 
^^C{i'e, e~y^X background, it further decreases the small background from the 
possible oscillation signal seen by LSND[^ 

Table || shows the efficiency for each of the selection criteria and the total 
efficiency for the process ^^C{ve, e~y^N* . Table shows the calculated number 
of background events from various sources satisfying these selection criteria. In 
addition to the background sources already discussed, there are several smaller 
sources of background. 

The neutral current excitation ^^C{v,iyy^C* (1+,1; 15.1 MeV) leads to 
prompt decay to photons with a 90% branching ratio. Most of the photons 
are eliminated by the 18 MeV energy requirement but approximately 0.8% will 
have reconstructed energies above 18 MeV due to the finite energy resolution. 
The measured cross section for this process Q is in good agreement with 
theoretical calculations!^ ^, ^. 

As discussed earlier, LSND has also measured the process ^^(7(1^^, /i^)^^iV* 
using a beam of j/^ from 7r+ DIF[^. Normally we detect both the /i~ and the 
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Table 6: The electron selection efficiency, the efficiency for satisfying the /3 
rejection criteria, the efficiencies for satisfying the angular and the energy re- 
quirements and the total efficiency for the process ^^C{i'e, e~)^^N* 



Quantity 


Efficiency 


Electron Selection 


0.497±0.035 


No Future /3 


0.997±0.001 


cos(? < 0.9 


0.980±0.005 


18 McV < E, < 34 UcY 


0.720±0.040 



Total Efficiency 0.349±0.031 



Table 7: Calculated number of background events satisfying the selection criteria 
for ^'^C{i'e,e~y^N* . Systematic uncertainties in the numbers of background 
events are also shown. 



Source 


Events 


'■'C{iy„e-y'Ng.,. 


434.1 ±28.1 


ue — > ve 


35.2 ±4.1 




46.6 ±23.3 


^'^Ciu,u)^'^C*[15.1 MeV] 


17.5 ±2.0 




9.3 ±2.0 


Events with a neutron 


21.5 ± 15.1 


Total Background 


564 ± 40 
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Table 8: Events, efficiency, neutrino flux and flux averaged cross section with 
statistical errors only for ^^C(i/e, e~)^^A^* 



Beam-on events 


2i26 


Beam-excess events 


1088 


Background 


564 




524 


Efficiency 


0.349 ±0.031 


Flux 


10.58 X 10^3 cni-2 




(4.3 ± 0.4) X 10-^^ cm^ 



Michel from the fi~ decay. Occasionally, however, we will miss a low energy 
li~ because it does not satisfy the activity threshold of > 17 PMT hits. In that 
case the Michel electron will be a background event if it satisfies the selection 
criteria given in Table ^. We have calculated this background using the event 
simulation program LSNDMC. We have also estimated this background using 
the observed distribution of PMT hits for /i^ satisfying the activity threshold. 

Finally, we determine the background from processes with associated neu- 
trons. Events from the reactions i?e + p + n and f'p + p — *■ fi'^ + n are 
identified in LSND by detecting the 2.2 MeV 7's from the capture reaction 
n + p ^ d + 7j|, |l3[ g3|. The mean capture time in the LSND detector is 186 
/xs, essentially independent of the initial neutron energy. Three variables are 
used to identify a capture 7 correlated with a neutron in the primary event: the 
number of PMT hits for the 7, the distance of the 7 from the primary event 
and the time of 7 from the primary event. A likelihood technique, discussed 
in Ref. has been developed to separate the correlated component due 
to neutrons from the uncorrelated or accidental component. An approximate 
likelihood ratio R — Ccor / C,uncor is calculated for each event from the three 
measured variables. We use here an analysis with improved photon spatial re- 
construction as is presented in Ref. ||3^, which reports final LSND results on 
neutrino oscillations. 

Figure |l^ shows the observed R distribution for events satisfying the selec- 
tion criteria for the process ^^C(i^e, e^Y'^N* . Shown for comparison is the best 
fit (solid line) combination of the correlated distribution and the uncorrelated 
distribution to the data. The best fit has a 2.0 ± 1.4% correlated component. 
The number of background events with neutrons shown in Table |^ is obtained 
from the best fit to the data. 

Table ^ shows the number of beam-on and beam-excess events satisfying the 
selection criteria, the number of background events and the resulting number 
of events and cross section for the process ^^C{ve,e.~Y^ N* . It is worth noting 
that the beam-off subtraction is much larger in this case than it was for the ex- 
clusive reaction where we require a coincidence with a (3. The flux average cross 
section is (cr) = (4.3 ± 0.4 ± 0.6) x 10"'*^ cm^. There are several contributions 
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Figure 15: The observed distribution of the 7 hkehhood ratio R for the process 
^'^C{i'e,e~y^N*. Shown for comparison is the best fit (sohd hne) combination 
of the correlated distribution and the uncorrelated distribution to data. The 
best fit has a 2.0 ± 1.4% correlated component. 
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Table 9: Measurements and theoretical predictions of the flux averaged cross 
section for the process ^^C(i/e, e~)^^iV*. 
Experiment 

LSND (4.3 ±0.4 ±0.6) x 10"^^ cm^ 

LSND(previous) § (5.7 ± 0.6 ± 0.6) x lO-'*^ 
E2250 111 (3.6 ± 2.0) x lO''^^ ^,^^^2 

KARMENH (5.1 ±0.6 ±0.5) x 10'"^^ cm^ 



Theory 



Kolbe eA al.\ 


% 


6.3 


X 


10" 


-42 


cm^ 


Kolbe eA al.\ 




5.5 


X 


10- 


-42 


cm^ 


Hayes et al.\ 




4.1 


X 


10- 


-42 


cm^ 



to the systematic error. The 7% flux uncertainty and the 6% uncertainty in the 
effective fiducial volume have been described previously. There is a 4% uncer- 
tainty due to the 50% error in the ^^C cross section. The uncertainty in the 6+ 
acceptance for the ^"^Ng^s. background subtraction leads to a 5% uncertainty in 
the ^"^N* cross section. The uncertainty in the duty ratio results in a 3% error 
in the cross section. We estimate a 5.5% uncertainty in the fraction of events 
with electrons in the region 18 MeV< Ee < 34 MeV. 

The cross section reported here is lower than that previously measured by 
LSND Hi primarily because of an increase in the calculated number of back- 
ground events, (see Table 0). An increase in the calculated deadtime due to 
the veto lowered the /3 selection efficiency and increased the calculated number 
of events from the process ^^C(j'e, e~)^^7Vg.s. without an identified f3. There 
was also an increase in the calculated number of events with cos 6* < 0.9 from 
the process lye — > i^e. The small DIE background from ^'^C{v^, fi^Y'^N* was 
not previously included. Finally, as was stated in Ref [|| , we previously did not 
subtract the background from events with neutrons. 

The flux averaged cross section measured by LSND is compared in Table || 
with other measurements and with several theoretical calculations. The most 
recent measurements of both LSND and KARMENp^ are somewhat lower than 
their previous measurements |j. The recent CRPA result of Kolbe et a/.[0, 
similarly, is lower than the previous result Q but it remains above the cross 
section measurements of both LSND and KARMEN. In contrast the recent 
shell model calculation of Hayes and Towner |^ is in good agreement with but 
lower than the cross section measured by LSND. 

Eigure |l6| shows that the measured electron energy distribution for the sam- 
ple of events satisfying the electron criteria given in Table || agrees well with 
that expected from ^^C(z/e, e~)^^iV* plus the background sources hsted in Ta- 
ble 0. Also shown is the energy distribution expected from just the background 
processes. Eor the largest background source, ^^C(i^e, e~)^^iVg.s., we used the 
shape of the energy distribution measured for the sample of events with identi- 
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Figure 16: The observed distribution of electron energy for events satisfying 
the selection criteria given in Table ^. The dotted (solid) histogram shows the 
expected distribution from all sources (background sources only). 
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Figure 17: The observed and expected (solid line) electron energy distribution 
for the process ^'^C{ve,e-Y'^N* 



fied /3's. For all other processes the expected energy distributions are obtained 
using simulated samples of events. Figure |l^ compares the measured energy dis- 
tribution after subtraction of the calculated backgrounds with the distribution 
expected for the process ^^C{Vf,,e~Y^N* . The agreement is excellent. 

Figure |l^ compares the measured distribution of cos6' after subtraction 
of the calculated backgrounds with the distribution expected for the process 
^'^C{Vf.^e~Y'^N* . The cos0 distribution is enhanced in the backward direction 
as expected. The data are in fair agreement with expectations j29| with the 
X^/DF = 1.79. The mean observed value of cos 6* is -0.15 ± 0.05(sta<.) ± 
0.030(syst.) compared with the expected value of —0.25 for Ref. |29) and —0.30 
for Ref. Q . The backward peaking of the angular distribution is largely a result 
of the negative parity of the N* states expected to contribute, 2~ levels at 1.20 
and 4.14 MeV and 1~ levels at 6.40 and 7.68 MeV. The I — 1 angular momen- 
tum transfer to the A — 12 system favors momentum transfer of approximately 
100 MeV/c, and hence the backward peaking. 

The total charged current cross section for interactions on ^^C is obtained 
by adding the cross sections measured for ^^C{i'e, e~)^^iVg.s. and '^'^C{i'e, e~)^'^N* 
The resulting flux averaged cross section for ^^C{i^e, s'Y^N is {a) = (13.2 ± 
0.5 ±1.3) X lO-''^ cm2. 
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8 Conclusions 



The process ^^C(i'e, e~)^^A^g s. has been measured with a clean sample of 733 
events. For this process the cross section calculations using empirical form 
factors are expected to be very reliable. The flux averaged cross section is 
measured to be (8.9 ±0.3 ±0.9) x 10"''^ cm^ in reasonable agreement with other 
experiments and theoretical expectations. The energy and angular distributions 
also agree well with theoretical expectations. 

The measurement of the process ^^C(i^e, e~)^^iV* is more diflicult, primarily 
due to the significant background subtraction required. The measured cross 
section of (4.3 ±0.4 ±0.6) x lO""''^ cm^ is in agreement with other measurements. 
It is in better agreement with the shell model calculation of Hayes and Towner 
[H than with the CRPA calculation of Kolbe et. al. , but is compatible with 
both models. 
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